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AN ANALYTICAL STUDY OF AIRPLANE -AUTOPILOT RESPONSE 
TO ATMOSPHERIC TURBULENCE 


By Waldo I. Oehman 
Langley Research Center 

SUMMARY 

An analytical study has shown that, with proper selection of feedback gains, an 
automatic control system can reduce excursions in altitude of a jet transport flying m 
turbulence without increasing structural loads. The control system uses feedback of 
attitude -angle and pitch-rate signals to the elevator to provide attitude hold and uses 
feedback of altitude and altitude-rate signals to the throttle to provide altitude hold. 
Without the altitude-hold control, altitude excursions of the airplane are sufficient, in 
mild turbulence, to cause the airplane to exceed its allowable altitude corridor. The 
altitude-hold mode is effective in reducing altitude excursions by about 50 percent. In 
extremely severe turbulence the airplane exceeds its limit normal-load factor with the 
control system in operation. 


INTRODUCTION 

The use of autopilots to reduce the response of airplanes to turbulence has been a 
controversial subject for many years. Earlier, it was thought that autopilots increase 
airplane response to turbulence and pilots were instructed to turn off the autopilot in 
rough air. Late in 1968, based on an extensive simulation study, it was recommended 
that an operational autopilot ("turbulence flight director") be developed for use during 
flight through turbulence (ref. 1). Recently, the use of autopilots during flight in turbu- 
lence has become more widely accepted. However, there have been only a few studies 
made by using various combinations of feedback quantities that give quantitative responses 
to both vertical and horizontal gusts in autopilot-controlled flight. Some exploratory 
results, obtained by using different analyses, are presented in references 2 and 3. 

The present study extends the study of reference 2 by consideration of more than 
one flight condition and by inclusion of filters and servos in the automatic control system. 
The objective is to determine what controls and what feedback quantities can be used to 
reduce excursions in speed and altitude from the desired trim condition without increas- 
ing structural loads on the airplane. 



Responses to atmospheric turbulence of a medium jet transport airplane are cal- 
culated for two altitudes and for several Mach numbers. The basic airplane and the air- 
plane operating with each of three variations of the automatic control system are studied. 
The calculated variables of motion are the root-mean-square (rms) values of airspeed, 
angle of attack, attitude angle, altitude, pitch rate, and normal-load factor. The responses 
are calculated for both vertical and horizontal components of turbulence. 


SYMBOLS 


[B(a;'3 matrix defined by equation (9) 

C D drag coefficient, 


C 


L 


lift coefficient, 

qS 


Cj^qS - Cq< sin aqS 
qS 


C L = ^ 

L qS 

C m pitching-moment coefficient, g^dung moment 

qSc 

C T thrust coefficient, - - - rust 

1 qS 

c wing mean aerodynamic chord, meters 


force along X-axis 
force along Z-axis 


f j(U),f 2 (U),f 3 (U) functions of gust velocity 

g gravitational constant, meters/second^ (lg = 9.807 m/sec^) 

M matrix of frequency -response functions 

h altitude, meters 

h c command altitude, meters 


K 0 ,Kq,Kt automatic-control gains (see fig. 2) 
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airplane radius of gyration, meters 
longitudinal scale of turbulence, meters 
Mach number 
moment about Y-axis 
mass of airplane, kilograms 

normal-load factor, Normal acceleration 

g 

free-stream dynamic pressure, newtons/meter^ 

wing area, meters2 

Laplace variable, 1/second 

time, seconds 

nondimensional time, — 

c 

velocity of local air mass (gust velocity), meters/second, U 

longitudinal component of U, meters/second 

airspeed, meters/second or knots 

lateral component of U, meters/second 

vertical component of U, meters/second 

angle of attack, radians 

change of air density with altitude, 1/meter, --^r 

P dh 

elevator deflection angle, radians 
command elevator deflection angle, radians 
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throttle deflection angle, radians 

attitude angle, radians 

command attitude angle, radians 

nondimensional Laplace variable, ^ s 

nondimensional mass, 

’ p Sc 

air density, kilograms/meter^ 

rms (root-mean-square) altitude, meters 

rms normal-load factor 

rms pitch rate, radians/second 

rms horizontal component of gust velocity, meters/second 

rms airspeed, meters/second or knots 

rms vertical component of gust velocity, meter s/second 

rms angle of attack, radians 

rms attitude angle, radians 

power spectral density of horizontal component of gust velocity, 
meters^/ second^ 

power spectral density of vertical component of gust velocity, 
meter s2/second^ 

spatial frequency, — , radians 
c 

circular frequency, radians/second 
reduced frequency, — , radians 
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Notation: 


' (Dot) first derivative with respect to time t 
“ (Double dot) second derivative with respect to time t 

’ (Prime) first derivative with respect to nondimensional time t’; for example, 

a* = — 
dt* 

" (Double prime) second derivative with respect to nondimensional time t' 

A perturbation quantity 

Derivatives are defined as follows: 


C D V 

8Cd 
= V U 

dV 

Cmo; : 

_ 8 c m 

da 

C Dq! 

_ 9C D 

da 

Cm a' 

_ 9 c m 

90! ' 

C Lqi 

_ 9C l 

da 

Cm 5 

°e 

_ 9 C m 

96 e 

C L a , 

_ 8£l 

da' 

c m 0 , 

_ 9 ^m 

V 

_ 8C l 
86 e 

C T y = 

dCrn 

: V - 

av 


8Cx 

C T 6t : 


C L e . 

__ xj 

de r 

' 8 fit 

C my 

_ v !Es 

9V 

C T p = 

9Crp 

p ~w 


THEORY 

Equations of Longitudinal Motion 

The derivation of the equations of longitudinal motion of a power-fixed airplane is 
well documented. (See ref. 4.) Since the equations used in the present study have terms 
involving the variation of thrust with altitude change and because the derivation of equa- 
tions with these terms is not well documented, a derivation is presented in the appendix. 
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The equations of airplane motion are linear differential equations in perturbation 
variables. The mass of a rigid airplane is assumed to be concentrated at the origin of a 
system of moving wind axes as depicted m figure 1. The X-axis is always alined with the 
velocity vector. Although the airplane may rotate with respect to the moving axes, it is 
assumed that changes in the moment of inertia Iy are negligible. Finally, the longitu- 
dinal motions are assumed to be uncoupled from lateral motions. In accord with these 
assumptions, the equations of longitudinal motion in Laplace notation are written as 
follows: 


^2/jX - C Tv cos a + C Dy ^f - (C L - Cq^JAcx + C L A8 + c^C T ^ + C T jcos a 


Ah 

c 




= -C T 6 t + fi(U) 
5 t 


C T y Sln « + 2 C l)^ + 


C x cos a + C La + (c La , + 2fjx Aa + fc L0t - 2/jJa A 8 
- c/3^C x ^ + Cq^sin a = -C^ ^e + ^2^ 


> W 


Cm v J Y :+ (Cm a + C mQ[1 \)Aa + (c m0t A - 2 X 2 |a 0 = -C m g^6 e + f 3 (U) 


Aa - A8 + = 0 


The desired equilibrium conditions are 


C x cos a = Cp 




= -C T sin a + 


L ~ T 


C m - 0 


8 = a 


PSV 2 } 


( 2 ) 


Equations (1) are made nondimensional in the following ways: 

(1) The force equations are divided by — SV^ 

A 

(2) The moment equation is divided by — ScV^ 

2 
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(3) Time derivatives are taken with respect to nondimensional time t' = — 

c 

The coefficients of the perturbation variables in equations (1) are evaluated at equilib- 
rium conditions. Control inputs are throttle deflection angle 6^ and elevator deflection 
angle 6 e . These inputs are given as transfer functions in the next section as equa- 
tions (3) and (4). The functions of gust velocity fj(U), f 2 (U), and f 3 (U) represent the 
turbulence input. 


Equations for Automatic Control System 

Two variations of the automatic control system are to be considered both separately 
and combined. One uses altitude and rate-of-climb signals to produce throttle changes. 
The purpose of this "autothrottle " is to maintain a desired altitude. The other uses 
attitude -angle and pitch-rate signals to produce elevator deflections. The purpose of this 
"pitch autopilot" is to maintain a desired attitude angle. A block diagram of the auto- 
matic control system is shown in figure 2 The pitch autopilot, with gains K g = 2.36 
and Kq = 3.2, is the attitude -hold mode of a conventional autopilot. 

The transfer function for the pitch autopilot, relating elevator deflection to attitude 
angle, is given in Laplace transform notation as follows: 



Band-pass filter Servo 


The gains Kq and Kq are adjustable. 

The transfer function for the autothrottle, relating change of thrust coefficient with 
altitude, is given in Laplace notation as follows: 


Cq 1 

— V— (a) 

Ah/c 



-0.004(^ X + 0.025) K t c 


4 — X + 11(4.4 — X + 1 




Time lag 


( 4 ) 


The gain Kf- is adjustable. This autothrottle is essentially the one used in reference 2. 


Analytical Representation of Atmospheric Turbulence 

Discussion of the random properties of atmospheric turbulence and the derivation 
of analytical representations are given in reference 5. The Von Karman one -dimensional 
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power-spectral-density representations are, for the horizontal component of gust velocity, 


*<■(<•>’) - °U 2 f 


1 + 1.339 


Lw' 


5/6 


and, for the vertical component of gust velocity, 

\2 

Ji.aaAi£L 

3 1 


^w( w ') - 0 w 2 


1 +f(l.339 W 


nc r- 


1 + 1.339 


Leo 


11/6 


(5) 


(6) 


The reduced frequency u>' is related to spatial frequency and circular frequency u> 
by 



c 


and 


U) - — co' 
c 

The longitudinal scale of turbulence L may be interpreted, physically, in various 
ways (see appendix B of ref. 5), and its value for atmospheric turbulence depends on alti- 
tude. One physical interpretation is that L is a rough measure of the longest distance 
by which two points in a turbulence field may be separated before the correlation function 
between the velocities becomes zero. The scale of turbulence is assumed to be 0.762 km 
for altitudes of 6.096 km and 10.668 km. Mean-square values of horizontal and vertical 
components of gust velocity are denoted by ct u ^ and a w 2 ( respectively. The power 
spectra <f> u (c<j') and <f> w (w') are plotted in figure 3. Use of statistical properties of 
atmospheric turbulence requires the application of random-process theory to obtain the 
statistical properties of the airplane response. 

Random -Process Theory 

The necessary mathematical development of random-process theory may be obtained 
from reference 6. An important result relates the power spectral densities of the input 
and the output of a linear system. The relationship may be stated as 

*o(w') - |H(u>-)| 2 <f-i(w') (7) 
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where |jl(co'f] is the frequency -response function of the linear system to a sinusoidal 
input. The function is the power spectral density of the input; and the function 

$ 0 (u>') is the power spectral density of the output. Application of equation (7) to prob- 
lems involving the response of airplanes to atmospheric turbulence is discussed thor- 
oughly in references 5, 7, and 8. 


For the problem considered, there is a frequency -response function for each of the 
variables AV/V, Aa, A 6, and Ah/c. These frequency-response functions are obtained 
from the equations of motion (see eqs. (1)) written in matrix form with the substitution of 
X. = ico'. The resulting matrix of coefficients, which contain the control transfer functions, 
is represented as the matrix [z(iu>'j]. The coefficients in the gust-velocity functions 
fj(U), f 2 (U), and f 3 (U) are the elements of the matrix (b(u)')]. Thus, 

AV 

V 



(j = 1, 2, 3, 4; k* 1, 2, 3, 4) 


( 8 ) 


Z 11 



- C<p cos a + 2/iico' 


z 12 

z 13 

z 14 

Z 21 

z 22 

z 23 

z 24 


" C D a " C L 
= C L 

= cjS^C-p^ + Crpjcos a - c lco’ 
= C-py sin a + 2Cl 


= C T cos « + C L a + ( C L a , + 2 ^J 1CJ ' 

= -cfi(Crj>p + Cq^sm a 
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Z 31 - c my 


z 32 _ Cm a + Cm ffl iW 

z 33 = w ’) + Cm 0' iw ’ + Cm 6 e a| 1w ’ 

z 34 = 0 

Z 41 = 0 

z 42 = 1 

z 43 = - 1 

and 

z 44 = uo' 

The complex matrix [B(u>'j] is, for umt vertical gust components, 

-(% - C l) 

C La + C T cos a + icu'C L q ^ 

‘(" C m a “ lw ’ c m a .) 

0 

and is, for umt horizontal gust components, 

C D V - C T V cos a 
sin a + 2Cl 

-C my 

0 




(9a) 


(9b) 
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The elements of the matrix [B(ct>’)] relate the changes in lift, drag, and pitching moment 
to changes in angle of attack (see eq. (9a)) and in airspeed (see eq. (9b)) caused by vertical 
and by horizontal gusts, respectively Vertical gusts cause a change in rms angle of 
attack given by a w /v. The change in rms airspeed caused by horizontal gusts is ct u /v. 

The required frequency-response functions are the elements of the matrix [H(a>')] 
m 

[h(u>')] = jz(aj')) 1 [b(u> , )J (10) 

Equation (10) is the solution of equation (8) for the variables AV/V, Aa, A0, and Ah/c. 

Substitution of equation (10) and either equation (5) or equation (6) into equation (7) 
yields the power spectral densities for each of the output variables AV/V, Aa, A0, and 
Ah/c. A combination of the frequency -response functions for Aa and A 6 gives the 
frequency-response function for airplane normal-load factor; and iw' times the 
frequency-response function for A 6 gives the frequency-response function for the pitch 
rate. 

Another result from random-process theory is 

r 00 

cfi = \ <f>(u>')du>' , (11) 

J 0 

which states that the mean-square value o 2 of a random variable is the integral of its 
power spectral density. 

Equations (5), (6), (8), and (11) are programed on a digital computer to obtain the 
rms values of the longitudinal-response variables of an airplane. These variables are 
a y> a a> a 6> a h> a q> anc * a n- Obviously, the integration in equation (11) can be car- 
ried out only to a fimte value of The upper limit for the integral is large enough to 
span the frequencies of all the airplane modes and large enough so that further increase 
does not appreciably increase the value of the integral. 

RESULTS OF CALCULATIONS 

Airplane data on mass and dimensions are listed in table I. The airplane is a 
medium jet transport about the size of a Boeing 720B. Pertinent aerodynamic data are 
listed in table II for altitudes of 6.096 km and 10.668 km and a Mach number range of 0.6 
to 0.85. Combinations of feedback gains used in the present study are as follows: 
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Pitch autopilot 

Autothrottle 

Combined pitch autopilot 
and autothrottle 


Kq 

K t 

K 0 

Kq 

K t 

0.59 

0.8 

0.0026 

0.59 

0.8 

0.0260 

1.18 

1.6 

.0260 




2.36 

3.2 

.1300 




4.72 

6.4 






Stability Calculations 

The stability of each airplane -autopilot combination is determined by examination 
of the sign of the real parts of the roots of the characteristic equation of each system. It 
is necessary to assure that all systems are stable because the analysis leading to equa- 
tion (7) is correct only if the system is stable. The airplane with the automatic controls 
not operating has an unstable phugoid for all flight conditions. This instability is caused 
by elimination of the major phugoid damping term Cj) from the equations of motion. 

(See the appendix.) A low-gain pitch-angle feedback (Kg = 0.236) is used to stabilize the 
phugoid. Most other combinations of airplane, automatic controls, and flight conditions 
are stable. At an altitude of 6.096 km and a Mach number of 0.85, one mode associated 
with the pitch autopilot (Kg = 4.72 and Kq = 6.4) is unstable. 

Responses to Turbulence 

Responses of the airplane-autopilot combinations to both umt rms horizontal and 
vertical components of turbulence (that is, cr w = a u = 1 m/sec) are calculated for com- 
binations of feedback gains and flight conditions. The responses to umt rms vertical gust 
components are presented for altitudes of 6.096 km and 10.668 km in figures 4 to 7; 
responses to umt horizontal components, in figures 8 to 11. Although the rms variations 
of the variables are plotted as functions of Mach number, it should be emphasized that the 
calculations are made at discrete Mach numbers with the assumption that the airplane is 
m equilibrium at those Mach numbers. The plots permit comparisons of the responses 
at the various Mach numbers (calculations are made at M = 0.6, 0.7, 0.8, and 0.85) and 
permit interpolation to obtain responses at intermediate Mach numbers. 

DISCUSSION 

The responses of the airplane to unit rms components of gust velocity (o w = <j u = 

1 m/sec) are useful for analysis of the automatic -control operation. Because the equa- 
tions are linear (see eqs. (5) to (7)), responses to any degree of turbulence may be 
obtained by use of appropriate multiples of umt o w and ct u (that is, of 1 m/sec). 
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Responses to Vertical Gusts 

Basic airplane .- The responses of the basic stabilized airplane (Kg = 0.236 and 
Kq = 0) to unit rms vertical gust components are presented in figure 4. The rms varia- 
tions of airspeed cry, angle of attack a a , attitude angle Oq, and pitch rate <7q are — 
relatively small for all flight conditions. However, the rms variations of altitude 
and normal-load factor cr n are considered to be excessively large This qualitative 
description of the responses of the basic airplane is sufficient for analysis of the per- 
formance of the automatic control system. The purpose of the analysis is to determine 
which feedback quantities can be used effectively to reduce excursions of the motion vari- 
ables from a desired trim condition without increasing structural loads on the airplane. 

Effect of pitch autopilot . - The responses of the airplane with the pitch autopilot 
operating (design gains Kg = 2.36 and Kq = 3 2) are also shown in figure 4. Since the 
purpose of the pitch autopilot is to maintain the attitude angle at some desired value, it is 
not surprising that the rms variation in attitude angle Oq is reduced as shown. This 
autopilot also produced sizable reductions of the rms variation of angle of attack o a , 
pitch rate Qq, and normal-load factor a n . The effect of the pitch autopilot on rms vari- 
ations of airspeed cry and altitude depends on Mach number and altitude. For any 
flight condition, the rms variation of altitude cr^ is not appreciably smaller than for the 
basic airplane; and, for most flight conditions, operation of the pitch autopilot results in 
an increased rms variation of airspeed Oy. 

Doubling the design gains of the pitch autopilot (that is, Kg = 4.72 and Kq = 6.4) 
does not produce an appreciable improvement in the response of the airplane to umt rms 
vertical gust components. (See fig. 5.) Proportionate decreases in the rms variations of 
attitude angle On and pitch rate Oq are obtained but the other variables, except air- 
speed, are only slightly increased or decreased. The rms variation of airspeed Oy gen- 
erally is proportionately increased. A proportionate increase or decrease of the rms 
value of a variable means that either doubling or halving the gains affects the value of the 
variable by approximately the same amount. 

When Kg = 4.72, Kq = 6.4, and h = 6.096 km, the rms variations of attitude angle 
Oq and pitch rate Oq have the unusual trend with Mach number shown by the curves 
labeled (4) in figure 5(a). In the section entitled "Stability Calculations," it is stated that 
a control mode is unstable at M = 0.85 and at h = 6.096 km. The damping of this mode 
decreases as the trim Mach number is increased. By way of explanation, the lightly 
damped mode can result in an increase of the power-spectral density at the natural fre- 
quency of the mode which, in turn, produces larger rms values of the variables. It should 
be recalled that the mean-square value of a variable is the integral of its power-spectral 
density. (See eq. (11).) 
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In practice, it is recommended that autopilot gains be reduced during flight in tur- 
bulence. Use of a pitch autopilot with larger than design gains cannot be recommended 
on the basis of the data presented herein. 

The rms variations of angle of attack o a , attitude angle Oq, pitch rate ct^, and 
normal-load factor a n increase when the pitch autopilot gains are reduced from the 
design values of K g = 4.72 and Kq = 6.4. However, the values of these rms variables, 
with Kg = 0.59 and Kq = 0.8, are smaller than the corresponding values for the basic 
airplane. (Compare curves labeled (1) of figs. 4 and 5.) Reducing the gains produces 
only small changes in the rms variation of altitude o^; but the rms variation of airspeed 
cry is reduced to approximately the values for the basic airplane. The pitch autopilot 
with gams Kg = 0.59 and Kq = 0.8 is considered to provide improved responses. 

Effect of autothrottle . - Responses to unit rms vertical gust components for the basic 
airplane and for the airplane with the autothrottle in operation are compared in figure 6. 
The rms variations of all the variables except airspeed <7y and altitude are prac- 
tically unchanged by operation of the autothrottle. Increasing the autothrottle gam K^ 
provides a large reduction in rms variation of altitude a^. The thrust changes required 
to hold altitude, however, cause an increase of the rms variation of airspeed Oy. The 
auto throttle gain Kj. = 0.0260 gives a sizable reduction of the rms variation of altitude 
cr h without a large increase of the rms variation of airspeed ay. So, K^ = 0.0260 is 
adjudged an improvement. 

Effect of pitch autopilot combined with autothrottle .- The responses of the airplane 
to vertical gusts with the combined pitch autopilot and autothrottle in operation (Kg = 0.59, 
Kq = 0.8, and K^ = 0.0260) are compared with the responses of the basic airplane in fig- 
ure 7. The rms variations of angle of attack a a , attitude angle <7g, pitch rate Oq, and 
normal-load factor a n with the combined controls in operation are essentially the same 
as for the operation of the pitch autopilot alone (Kg = 0.59 and Kq = 0.8). (See fig. 5.) 
The rms variations of airspeed ay and altitude with the combined controls in oper- 
ation are essentially the same as for the operation of the autothrottle alone (Kt = 0.0260). 
(See fig, 6.) However, the rms variation of airspeed cry at an altitude of 6.096 km is 
much larger with the combined controls m operation than with the autothrottle alone in 
operation. 


Responses to Horizontal Gusts 

The responses of the airplane to unit rms horizontal gust components (ct u = 1 m/sec) 
are characterized by relatively small rms variations of the variables a a , Oq, o^, Oq, 
and o n . (See figs. 8 to 11.) The rms variation of airspeed Oy is as large or larger 
than for the response to umt rms vertical gust components. The rms variation of altitude 
is less than 10 meters for all flight conditions but, in severe turbulence, this variation 
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may be more than desired. In figures 8 to 11, the rms variations of pitch rate crq are 
not accurately represented. The calculated values of are so small that, without a 
change of scale, the curves lie almost entirely on the axis at Oq = 0. Therefore, the 
lines shown m figures 8 to 11 represent the largest calculated value of Oq. The rms 
variations of angle of attack o a , attitude angle Oq, and pitch rate Oq are almost neg- 
ligible even m severest turbulence. The effect of the automatic control system on the 
responses to unit rms horizontal gust components is similar to the effect of the automatic 
control system on the responses to unit rms vertical gust components. 

Assessment of Control Systems 

In the responses of the airplane to unit rms vertical gust components, the automatic 
control system with selected feedback gains produces the desired reduction in excursions 
of altitude without increasing the normal-load factor. However, excursions m airspeed 
are increased. The pitch autopilot is not effective in reducing excursions in altitude; and 
the autothrottle, although quite effective in reducing excursions in altitude, causes 
increased excursions in airspeed. 

The relevance of the control effectiveness is shown by consideration of a practical 
situation. It is assumed that the airplane is trimmed for M = 0.85 at h = 10.668 km 
and is flying in severe turbulence (that is, cr w = 9 m/sec). The proximity to the buffet 
boundary for the basic airplane is illustrated in figure 12(a) by the parallelogram sur- 
rounding the trim condition. The rms variation of altitude is = 297 m and the rms 
variation of airspeed is <7y = 1.7 knots The curve labeled n = 1 53 is the buffet bound- 
ary for a normal-load factor of 1.53. (This corresponds to the calculated rms variation 
of normal-load factor a n = 0.53.) With the pitch autopilot in operation (kg = 0.59 and 
Kq = 0.8), the airplane response to severe turbulence is illustrated by the parallelogram 
with solid lines in figure 12(b). The rms variation in altitude is increased to ojj = 333 m 
and the rms variation of airspeed is increased to cry = 1.8 knots. However, the rms 
normal-load factor is reduced to a n = 0.46. The effect of adding the autothrottle = 
0.0260) is illustrated by the parallelogram with dashed lines m figure 12(b). The rms 
variation of normal-load factor remains the same; but the rms variation of altitude is 
reduced to 0 ^ = 144 m or about one-half that of the basic airplane and the rms airspeed 
is further increased to <7y = 3.4 knots. 

In this severe turbulence (a w = 9 m/sec), <jy, a h , and a n remain below the cor- 
responding buffet boundary for the airplane, even with the automatic controls operating. 
Extreme gust velocities, however, can be three times their rms values Thus, with alti- 
tude excursions between about 430 m and 1000 m, depending on the control system, the 
buffet boundary is crossed. Also, the limit normal-load factor of 2.5 is exceeded. 
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The responses of the airplane shown in figure 12 for a w = 9 m/sec may be con- 
sidered to represent the extreme condition, for which the gust component is three times 
its rms value of 3 m/sec. The value of ct w = 3 m/sec characterizes mild turbulence. 

It is not likely that the buffet boundary is crossed while the airplane is flying in this mild 
turbulence; but the altitude excursions, without altitude hold, are sufficient to cause the 
airplane to exceed its allowable altitude corridor of 600 m. Thus, the danger of collision 
is increased. Further, the wide variations of normal-load factor can be uncomfortable 
for passengers. 


CONCLUDING REMARKS 

The present study has shown that, with proper selection of feedback gains, an auto- 
matic control system can reduce excursions in altitude of a jet transport flying in turbu- 
lence without increasing structural loads. The control system uses feedback of attitude- 
angle and pitch-rate signals to the elevator to provide attitude hold and uses feedback of 
altitude and altitude-rate signals to the throttle to provide altitude hold. Without the 
altitude-hold control, altitude excursions of the airplane are sufficient, in mild turbulence, 
to cause the airplane to exceed its allowable altitude corridor. The altitude-hold mode is 
effective m reducing altitude excursions by about 50 percent. In extremely severe turbu- 
lence the airplane exceeds its limit normal-load factor with the control system in 
operation. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., July 26, 1972. 
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APPENDIX 


EQUATIONS FOR AIRPLANE LONGITUDINAL MOTION 

It is assumed that the mass, a constant, of an airplane is concentrated at the origin 
of a system of moving wind axes as shown in figure 1. Changes in the moment of inertia 
about the Y-axis of the airplane Iy = mky 2 are assumed to be negligible. Aerodynamic 
and applied force components act along the X- and Z-axes with a moment about the Y-axis 
The airplane is free to move in the X- and Z -directions and to rotate about the Y-axis. 
The equations of longitudinal motion, starting from equilibrium, are obtained by expand- 
ing, in a Taylor's series, the following equations describing Newton's laws of motion: 

f x (V,V,a,0,h,U,6 t ) = F X - mV = 0 

f Z (V,V,a,0,h,a,0,6 e ,U,6t) = F z + mV (0 - d) = 0 

M y (V,V,a,d,0,U,6 e ) = My - mky 2 0 = 0 


Dependence of the variables on time is implicit in these equations. In the series expan- 
sion, the first terms and the coefficients of the remaining terms are evaluated at equilib- 
rium conditions. Also, forces and moments caused by control movement and by atmo- 
spheric turbulence are placed on the right-hand side of the equations. An equation for 
rate of climb is given since force variations with changes in altitude are included. * For 
small disturbances, only the first two terms of the series need be retained. The result- 
ing linear equations are given as follows: 


• 9F X av 9F X 9F X 9 F X Ah AV 

F x - mV + V -# ^ Aa + — £ A 0 + c — ^ ^ - mV = -F x ( 6 t ,U) 

A aV V 0a 00 ah c V A 1 


9F Z AV 9F Z 


9F Z A a - 9F Z Ah 


F Z + mV(0 - d) + V — f ^ Aa + — P- A a + — - A 0 + — ^ A 0 + c — - ^ 

41 0V V 0a 0a 00 90 0h c 


+ mV(0 - a)^ + mV(A 0 - Ad) = -F z ( 6 e ,U) 


9 My 0 My 


M v - mkv^ + V — - — + — - A a + — ^ Ad + — - A 0 - mk Y 2 A 0 '= -My( 6 e ,U) 
Y Y 0V V 0a 0d 80 Y 1 


h - V sin (0 - a) - V sin (0 - a) — — + V cos (0 - a)Aa - V cos (0 - a)A 0 + c — = 0 

V c 
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APPENDIX - Continued 


Equilibrium conditions are characterized by zero acceleration and zero rate of 
climb (that is, V = V0 = Va = 0 = h = 0). This condition leads to 

F x = 0 
F z = 0 
My = 0 


6 = a 


Also, it is assumed that F^(6t,U) = F z (6 e ,U) = My(6 e ,U) =0 at equilibrium. 
Explicitly, the forces and moments are given by the following equations: 


F x = - SV 2 


F z = § SV 2 


jc T (V,h)cos a - Cj)(o!,V) - C L sin(0 - ajj 
|jc-j’(V,h)sin a - CL(a,a,fl) + Cl cos(0 - a) 


3 


My = — ScV 2 C m (V,Qf ,d,0) 


J 


and the equilibrium conditions are as follows: 

Cj) = Cy COS OL 

C L = C L - C T sin a \ 


(A2) 


(A3) 


The desired linear differential equations are obtained by use of equations (A3) and 
by substitution of equations (A2) into equations (Al). The desired forms of the equations 
include the following: 

(1) After substitution, force equations are divided by £ SV 2 and moment equations, 

z 

by |S5V 2 

(2) Time t is changed to nondimensional time t' = ^ with derivatives with 

c 

respect to t' denoted by a prime and thus, for instance, — = — — — = — a' 

dt 5 tV c 

c 
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APPENDIX - Concluded 


(3) Mass is nondimensionalized as p = 

(4) Change of air density with altitude is j3 = which comes from p = P Q e~^ 

when p Q is air density at standard sea-level conditions 


The resulting equations are as follows: 


2 M - (c Ty cos a - + 

+ C-j^cos a — = -C T 6t + f l( u ) 


( C D a ~ 


CjjAa + Ct. A 6 + BciCr 


( C T V sin a + 2 C l)^ + ( c L a + C T cos ot)&a + ^2p + C^Acf + ( C L0 , - 2p)&6' 
- c/S^C-pp + C-pjsin a = -CL^dg + f2(U) 


Cmy ~y~ + Cm a ^ a + C mQ ,t Aa' + Cmgi A0' “ 2p — — A0" = "^mg^e + f3(U) 


Aa - A0 + — = 0 


For power-fixed flight, the coefficient of AV/V in the first of equations (A4) is 
+(C D + Cj)^. It is this term that provides most of the damping of the phugoid. The 
Cp-part of the term is eliminated by the thrust component for powered flight and, conse- 
quently, the phugoid is likely to be unstable without some means of stability augmentation. 
The Laplace transforms of equations (A4) are taken with the Laplace variable A. = ^ s to 
obtain the following equations: 
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TABLE I.- AIRPLANE DATA 


Mass, kg 79 379 

Wing area, 236.9 

Mean aerodynamic chord, m 5.86 

Wing span, m 40.69 

Moment of inertia about Y-axis, kg-m^ 4.74 
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TABLE II.- AIRPLANE AERODYNAMIC DATA 
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Figure 2.- Block diagram of automatic control system. 




















(a) Altitude = 6.096 km. 

Figure 5.- Effect of pitch-autopilot gain on responses of airplane to vertical gusts. a w = 1 m/sec. 
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Mach no Mach no Mach no 



(a) Altitude = 6.096 km. 

Figure 6.- Effect of autothrottle on responses of airplane to vertical gusts. a w = 1 m/sec. 



(I)(2)(3) (4) 



Figure 6.- Concluded. 





Figure 7.- Effect of combined pitch autopilot and autothrottle on responses of airplane to 

vertical gusts. cr w = 1 m/sec. 



Figure 7.- Concluded. 








Figure 8.- Concluded. 






Figure 9.- Effect of pitch-autopilot gain on responses of airplane to horizontal gusts o u = 1 m/sec. 




Figure 9.- Concluded. 






(a) Altitude = 6.096 km. 

Figure 10.- Effect of autothrottle on responses of airplane to horizontal gusts. a„ = 1 m/sec. 





Figure 10.- Concluded. 








Figure 11.- Concluded. 
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Indicated airspeed, knots Indicated airspeed, knots 

(a) Basic airplane. Kg = 0.236; Kq = 0.0. (b) Airplane with automatic controls. Kq = 0.59; 

Kq = 0.8; K t = 0.0260. 

Figure 12.- Proximity to buffet boundary during flight in turbulence (vertical gusts), a = 9 m/sec. 
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